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Topological insulators are materials with a bulk excitation gap generated by the spin orbit inter- 
action, and which are different from conventional insulators. This distinction is characterized by 
Z2 topological invariants, which characterize the groundstate. In two dimensions there is a single 
Z2 invariant which distinguishes the ordinary insulator from the quantum spin Hall phase. In three 
dimensions there are four Z2 invariants, which distinguish the ordinary insulator from "weak" and 
"strong" topological insulators. These phases are characterized by the presence of gapless surface 
(or edge) states. In the 2D quantum spin Hall phase and the 3D strong topological insulator these 
states are robust and are insensitive to weak disorder and interactions. In this paper we show that 
the presence of inversion symmetry greatly simplifies the problem of evaluating the Z2 invariants. 
We show that the invariants can be determined from the knowledge of the parity of the occupied 
Bloch wavefunctions at the time reversal invariant points in the Brillouin zone. Using this ap- 
proach, we predict a number of specific materials are strong topological insulators, including the 
semiconducting alloy Bii-^ Sb x as well as a-Sn and HgTe under uniaxial strain. This paper also 
includes an expanded discussion of our formulation of the topological insulators in both two and 
three dimensions, as well as implications for experiments. 

PACS numbers: 73.43.-f, 72.25. Hg, 73.20.-r, 85.75.-d 
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I. INTRODUCTION 

In elementary solid state physics textbooks an insula- 
tor is described as a material with an energy gap separat- 
ing filled and empty energy bands. A more sophisticated 
definition of an insulator is that of a material for which 
all electronic phenomena are loca^. Such a material is 
insensitive to boundary conditions, so that in a multiply 
connected sample, such as a ring, there is exponentially 
small sensitivity to magnetic flux threading the holes. 
The simplest atomic insulator, in which all electrons are 
tightly bound to atoms in closed shells, clearly satisfies 
both properties. Ionic and covalent insulators do too. 
These band insulators are topologically equivalent in the 
sense that the Hamiltonian can be adiabatically trans- 
formed into an atomic insulator without going through 
any phase transitions. With regards to their low energy 
electronic behavior, conventional insulators are equiva- 
lent to atomic insulators. 

The existence of a bulk energy gap does not guarantee 
the insensitivity to boundary conditions, and there exist 
phases with bulk gaps, which are topologically distinct. 
In addition to exotic strongly correlated phases^ 3 -, this 
fact arises even for non interacting electrons described 
within band theory. The simplest example is the inte- 
ger quantum Hall effect (IQHE). In two dimensions, a 
magnetic field introduces a cyclotron gap between Lan- 
dau levels, which may be viewed as energy bands in the 
magnetic Brillouin zone. This phase can exist even with- 
out Landau levels in the absence of a uniform magnetic 
field 4 -, though a necessary condition is that time reversal 
symmetry be broken. Based on the bandstructure alone 
it is difficult to tell the difference between the IQHE state 
and a band insulator. The distinction between the two 
is a topological property of the occupied bands which is 
encoded into the Chern integer introduced by Thouless 



et al£. Three dimensional generalizations of the IQHE 
state, which may be viewed as layered 2D states, are in- 
dexed by a triad of Chern integers^. A hallmark of the 
IQHE phases, which is intimately related to their topol- 
ogy, is the existence of gapless chiral edge states which 
are robust in the presence of disorder^. 

Recently, new topological insulating phases for 
systems with time reversal symmetry have been 
discovere d 9 ! 10 ' 11 ! 12 ' 13 ! 14 ' 15 . In two dimensions, the quan- 
tum spin Hall phase is distinguished from a band insula- 
tor by a single Z2 invariant^. This phase exhibits gap- 
less spin-filtered edge states, which allow for dissipation- 
less transport of charge and spin at zero temperature, 
and are protected from weak disorder and interactions 
by time reversal symmetry. In three dimensions a time 
reversal invariant bandstructure is characterized by four 
Z2 invariant s 13 ' 14 ' 15 . Three of the invariants rely on the 
translational symmetry of the lattice and are not robust 
in the presence of disorder, leading to "weak topological 
insulators" . The fourth invariant, however, is robust and 
distinguishes the "strong topological insulator" (STI). 

Nontrivial Z2 invariants imply the existence of gapless 
surface states. In particular, in the STI phase, the sur- 
face states form a two dimensional "topological metal" , 
in which the Fermi arc encloses an odd number of Dirac 
points^ 5 .. This leads to a quantized Berry's phase of w ac- 
quired by an electron circling the surface Fermi arc, which 
does not change under continuous perturbation s 16 ' 17 . 
The 7r Berry's phase also signifies that with disorder the 
surface states are in the symplectic universality class, 
and exhibit antilocalization— . Thus, the metallic surface 
states form a unique phase, which can not be realized in 
a conventional two dimensional electron system for which 
Dirac points must come in pairs^ 9 -. 

The purpose of this paper is twofold. First, we will 
explain the formulation of the Z2 invariants in somewhat 
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more detail than in Ref. [HI . emphasizing the physical 
meaning of the invariants and their relation to the sur- 
face states. Second, we will develop a new method for 
evaluating the Zi invariants in systems which have in- 
version symmetry. This allows us to predict a number of 
specific materials to be strong topological insulators. 

Most insulators are conventional insulators. In order 
to find topological insulators experimentally it is neces- 
sary to develop criteria for recognizing them from their 
bulk band structure. Clearly, a necessary condition is the 
existence of a bulk energy gap which owes its existence 
to the spin-orbit interaction. However, evaluating the Zi 
invariants for a given band structure is in general a dif- 
ficult problem. To date three general approaches have 
been used. 

(1) One powerful approach is to exploit adiabatic con- 
tinuity to a Hamiltonian which has extra symmetry. We 
used this method to identify the quantum spin Hall 
phase in graphene^^ by arguing that the Hamiltonian 
for graphene is adiabatically connected to a Hamiltonian 
in which the spin S z is conserved. With this extra con- 
servation law the system can be characterized by a spin 
Chern number, which describes the quantized spin Hall 
conductivit y 20 ' 21 . The Z2 invariant can then be identified 
with the parity of the spin Chern number. In the presence 
of S z non conserving terms (which are inevitably present) 
the spin Chern number loses its meaning 22 . However, the 
Z2 invariant retains its value and characterizes the quan- 
tum spin Hall phase. 

Adiabatic continuity can also be used to establish a 
material is a band insulator if an adiabatic path can be 
found which connects the material to an "atomic" limit. 
Moreover, it can be argued that the Z2 invariant changes 
at an appropriate quantum phase transition, where the 
bulk energy gap goes to zero 12 ' 14 . In general, this ap- 
proach require a continuous path be found which con- 
nects the Hamiltonian in question to a known phase. 

(2) It is also possible to evaluate the Z2 invariant di- 
rectly with the knowledge of the Bloch wavefunctions for 
the occupied energy bands. In Ref. [H we established a 
formula for the invariant in terms of a Brillouin zone in- 
tegral. This is analogous to the calculation of the Chern 
number as an integral of the gauge invariant Berry's 
curvature^^. However, unlike the Chern invariant, the 
integral for the Z2 invariant also involves the Berry's po- 
tential, and requires a gauge in which the wavefunctions 
are globally continuous. Since time reversal symmetry 
requires the Chern invariant to vanish, a globally con- 
tinuous gauge is guaranteed to exist. However, finding a 
continuous gauge is not always simple. 

(3) A third approach is to characterize the zeros of 
Pfaffian function introduced Ref. Though the Pfaf- 
fian is not gauge invariant, its zeros can be determined 
without specifying a continuous gauge. While this ap- 
proach is tedious (especially in three dimensions) it has 
been successfully implemented by Murakami 2 ^ to show 
that 2 dimensional bismuth bilayers realize a quantum 
spin Hall phase. 



In this paper we will show that the presence of inver- 
sion symmetry greatly simplifies the problem of identify- 
ing the Z2 invariants. We show that the invariants can 
be determined from the knowledge of the parity of the 
occupied band eigenstates at the eight (or four in two 
dimensions) time reversal invariant momenta Ti in the 
Brillouin zone. Specifically, we will show that the Z2 
invariants are determined by the quantities 

JV 

Si = IJ &m(ri). (1.1) 

m— 1 

Here ^2m(r«) = ±1 is the parity eigenvalue of the 2m' th 
occupied energy band at Ti , which shares the same eigen- 
value ^2m — £,2m-i with its Kramers degenerate partner. 
The product involves the 2N occupied bands. The Z 2 
invariant v = 0, 1 which distinguishes the quantum spin 
Hall phase in two dimensions and the strong topological 
insulator in three dimensions is then given by the product 
of all the 5i's, 

(-ir=n^- (!-2) 

i 

The other three "weak" topological invariants in three 
dimensions also are determined by Si. Since the par- 
ity eigenvalues £ n (Ti) are tabulated in the band theory 
literature this allows us to identify inversion symmet- 
ric topological insulating materials. Moreover, exploit- 
ing adiabatic continuity allows us to identify topological 
insulators which don't have inversion symmetry, but are 
adiabatically connected to materials which have inversion 
symmetry. 

Applying the above approach, we predict that the fol- 
lowing narrow gap semiconductors are strong topological 
insulators: (1) the alloy Bii-^Sb^, which is semiconduct- 
ing for .07 < x < .22. (2) a— Sn and HgTe under uniaxial 
strain, (3) The alloy Pbi-^Sn^Te under uniaxial strain 
for x ~ x c in the vicinity of the band inversion transition. 
The materials (2-3) were suggested by Murakami, Na- 
gaosa and Zhang 2 ^ as candidates for spin Hall insulators. 
Those authors argued that those materials share a large 
spin-Hall conductivity, as calculated by a Kubo formula. 
Our analysis of these materials is rather different, and we 
will show that PbTe is a conventional insulator, despite 
its large spin-Hall conductivity, while strained a— Sn and 
HgTe are topological insulators. 

In section II we will present an expanded discussion of 
our formulation of the Z2 invariants. Then, in section III, 
we will derive Eqs. (1.1,1.2) for problems with inversion 
symmetry. In section IV we will apply our method to a 
class of four band tight binding models, which includes 
the graphene model as well as the 3D model introduced 
in Ref. 15. In section V we will apply (1.1,1.2) to deduce 
the Z2 invariants of several real materials based on their 
known band structures. Readers uninterested in the tech- 
nical details can skip directly to section V to read about 
these applications. Finally, in section VI we will conclude 
with a brief discussion of the experimental implications 
for the topological insulating phases. 
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II. Z-2 INVARIANTS IN TWO AND THREE 
DIMENSIONS 

In this section, we will review our formulation of the 
topological insulating phases. We begin in II A by defin- 
ing the time reversal polarization. In IIB, we develop the 
Z2 characterization of a bandstructure as a topological 
property of the occupied Bloch wavefunctions. In IIC we 
show how the i? 2 invariants determine the surface state 
spectrum. In IIC, we consider a more general formula- 
tion of the Zi invariant as a sensitivity of a bulk crystal 
to boundary conditions. 



A. The Time Reversal Polarization 

In Ref. [22] we introduced the concept of the time rever- 
sal 'polarization, in the same spirit as charge polarization. 
This allows for an interpretation of the Z2 invariants us- 
ing a Laughlin type gedanken experiment on a cylinder 26 . 
In the ordinary quantum Hall effect, the Chern invariant 
can be interpreted in a similar way. To motivate the time 
reversal polarization we therefore begin with a discussion 
of the charge polarization. 

The charge polarization determines the surface charges 
present in a finite system. Since electrons may be added 
or removed from a surface, the charge polarization is only 
defined modulo an intege r 27 i 28 i 29 ' 30 . However, changes in 
the charge polarization induced by adiabatic changes in 
the Hamiltonian arc well defined. In Laughlin's gedanken 
experiment for the integer quantum Hall effect a quan- 
tum of magnetic flux h/e is adiabatically inserted in a 
cylindrical quantum Hall sample at filling v = N. The re- 
sulting transfer of N electrons from one end of the cylin- 
der to the other can be interpreted as a change in the 
charge polarization of the cylinder. The Chern invari- 
ant, which characterizes the integer quantum Hall state 
precisely characterizes this quantized change in charge 
polarization. 

The time reversal polarization is a Z2 quantity, which 
characterizes the presence or absence of a Kramers de- 
generacy associated with a surface. Like the charge po- 
larization, its value can be changed by adding an extra 
electron to the surface. Thus by itself, the time rever- 
sal polarization is not meaningful. However, changes in 
the time reversal polarization due to adiabatic changes 
in the bulk Hamiltonian are well defined. Specifically the 
change in the time reversal polarization when half a flux 
quantum h/2e is threaded through a cylinder defines a 
Z2 invariant which is analogous to the Chern invariant, 
and distinguish topological insulators. 



B. Z2 invariants characterizing a bandstructure 

The topological invariant characterizing a two dimen- 
sional band structure may be constructed by imagining 



a long cylinder whose axis is parallel to a reciprocal lat- 
tice vector G, and which has a circumference of a single 
lattice constant. Then the magnetic flux threading the 
cylinder plays the role of the circumferential (or "edge" ) 
crystal momentum k x , with $ = and $ = h/2e corre- 
sponding to two "edge" time reversal invariant momenta 
k x = Ai and k x — A 2 . The Z2 invariant characterizes 
the change in the Kramers degeneracy at the ends of this 
I dimensional system between k x = K\ and k x = A 2 . 

For a three dimensional crystal, imagine a "generalized 
cylinder" which is long in one direction (parallel to G), 
but in the other two directions has a width of one lattice 
constant with periodic boundary conditions. While this 
structure can not be pictured as easily as a regular cylin- 
der, a distorted (but topologically equivalent) version can 
be visualized as a torus with a finite thickness. This 
"corbino donut" is analogous to the generalized cylin- 
der in the same way the corbino disk is analogous to the 
regular cylinder. The "long" direction corresponds to the 
thickness of the torus, and the two ends correspond to the 
inner and outer surfaces. This system can be threaded by 
two independent magnetic fluxes, and they correspond to 
the two components of the momentum perpendicular to 
G. There are four time reversal invariant surface mo- 
menta A a , corresponding to the two fluxes being either 
or h/2e. The bandstructure can be characterized by the 
difference in the time reversal polarization between any 
pair. 

The Z2 invariants can be deduced from the topological 
structure of the Bloch wavefunctions of the bulk crystal 
in the Brillouin zone. Consider a time reversal invariant 
periodic Hamiltonian H with 2N occupied bands charac- 
terized by Bloch wavefunctions 

|V>n,k) = e lk ' r Kk). (2.1) 

Here \u n ^) are cell periodic eigenstates of the Bloch 
Hamiltonian, 

H{k) = e~ lkr He +lkr . (2.2) 

We require lipn^+c) = IV'n.k) for reciprocal lattice vec- 
tors G, so that the Brillouin zone in which k is defined is a 
torus. This implies |u n ,k+G) = e~ lG ' r \u n ^). Time rever- 
sal symmetry implies [H, 6] = 0, where O = exjp(iTrS y )K 
is the time reversal operator (S y is spin and K com- 
plex conjugation), which for spin 1/2 particles satisfies 

e 2 = -1. it follows that H(-k) = eH(k)e-\ 

There are special points k = T,; in the Brillouin zone 
which are time reversal invariant and satisfy — Tj 
I , + G for a reciprocal lattice vector G. There are eight 
such points in three dimensions and four in two dimen- 
sions. At these points H(Ti) = 0_ff(r;)0 _1 , so that the 
eigenstates are Kramers degenerate. In the following it 
will be useful to use two different notations to label the 
distinct IV (1) The eight (or four) I\ can be indexed by 
three (or two) integers nj = 0, 1 defined modulo 2, which 
specify half a "mod 2 reciprocal lattice vector" , 

r i= ( ni „ 2 „ 3 ) = - (nibi + n 2 b 2 + n 3 b 3 ) , (2.3) 
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FIG. 1: (a) A two dimensional cylinder threaded by magnetic 
flux When the cylinder has a circumference of a single 
lattice constant $ plays the role of the edge crystal momen- 
tum k x in band theory, (b) The time reversal invariant fluxes 
$ = and h/2e correspond to edge time reversal invariant 
momenta Ai and Aa- A a are projections of pairs of the four 
bulk time reversal momenta r i= ( aft ), which reside in the two 
dimensional Brillouin zone indicated by the shaded region, (c) 
In 3D the generalized cylinder can be visualized as a "corbino 
donut", with two fluxes, which correspond to the two com- 
ponents of the surface crystal momentum, (d) The four time 
reversal invariant fluxes $i, <E>2 = 0, h/2e correspond to the 
four two dimensional surface momenta A a . These are projec- 
tions of pairs of the eight I\ = ( aM ) that reside in the bulk 3D 
Brillouin zone. 



where b; are primitive reciprocal lattice vectors. Two 
mod 2 reciprocal lattice vectors are equivalent if they 
differ by twice a reciprocal lattice vector. (2) A second 
notation is useful when considering a surface perpendic- 
ular to reciprocal lattice vector G. The surface time re- 
versal invariant momenta A Q will be projections of pairs 
of Ti which differ by G/2, as shown in Fig. [TJ Given G, 



we can define T, 



= (a/i) 



such that r a 



G/2. 



The change in the time reversal polarization associated 
with a cylinder oriented along G due to changing the 
flux between two values corresponding to A a and A& can 
be calculated by a method analogous to the calculation 
of the c harg e polarization as a Berry's phas o 27 ' 28 i 29 ' 30 . 
In Ref. [23 we showed that the result is most simply 
expressed in terms of the quantities 



Si = 



= ±1, 



where w is the 2A*" x 27V unitary matrix defined by 
uw(k) ee (u m _ k |9|u„k)- 



(2.4) 



(2.5) 



Since (Qa\Qb) = (b\a) and 9 2 = —1 the matrix w(k) 
is antisymmetric for k = Tj, allowing for the definition 
of the Pfaffian, which satisfies detky] = Pf[w] 2 . Given a 
surface G the time reversal polarization associated with 
surface momentum A„ was found to be 



Sal$a2 



(2.6) 



This formula is analogous to the expression for the charge 
polarization as a Berry's phase along a closed cycle in 
momentum spaced. 

Because of the square root, the sign of Si is ambigu- 
ous. However, since we require |u rai k) to be continu- 
ous, yj det [w(k)] is defined globally throughout the Bril- 
louin zone. The sign ambiguity then cancels for 7r a . 
The existence of the global square root is guaranteed 
by the fact that the phase winding of det[u;(k)] around 
any closed loop C is zero, so that -y/det[w(k)] has no 
branch cuts. When C is contractable, the vanishing 
phase winding follows from the continuity of \u n ^}. For 
non contractable loops, which can be continuously de- 
formed to satisfy C = —C, it follows from the fact that 
det[iu(-k)] = det[iu(k)]. 

The inherent ambiguity of the time reversal polariza- 
tion is reflected in the fact that the 7r a are not gauge 
invariant. For instance, consider a gauge transformation 
of the form 



|Wn,k) 



e i9k \u, 



= 1 



,k) for n 
, )k ) for n^l. 



Under this transformation, 



det[iu(k)] 
Pf[«j(r,)] 



det[u)(k)]e" 
Pi[w(Ti)]e 



(2.7) 



(2.8) 
(2.9) 



If we choose #k = k-R for a lattice vector R the Brillouin 
zone periodicity of u„k is preserved. From [2~51 it is clear 
that det[u>(k)] is unchanged. However, if G • R = 27r, it 
follows that Pf [u>(r o i)]Pf [u>(r a 2)], and hence 7r a , changes 
sign. Since this gauge transformation changes the sign 
of 7r Q for all a, however, the product 7r a 7rt,, which gives 
the change in the time reversal polarization between A a 
and A?, remains invariant. In general, the product of any 
four 5j's for which lie in the same plane is gauge in- 
variant, and defines a topological invariant characterizing 
the band structure. 

In two dimensions there are four time reversal invariant 
momenta Ti and a single Zi invariant, given by 



(-ir = n* 



(2.10) 



In three dimensions there are 8 time reversal invariant 
points. This leads to 4 independent Z-i invariant o 13 i 14 ' 15 . 
One of these invariants, vq, can be expressed as the prod- 
uct over all eight points, 



(-ir=n^ 



(2.11) 



The other three invariants are given by products of four 
Si's for which reside in the same plane. 



(-ir = n 



^ i—(riiri2ri^) • 



(2.12) 
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vq is clearly independent of the choice of primitive re- 
ciprocal lattice vectors bj, in (|2.3[) . {y\ViV%) are not. 
However, they may be viewed as components of a mod 2 
reciprocal lattice vector, 



(2.13) 



This vector may be explicitly constructed from the Si's 
as follows. A gauge transformation of the form (|2.7[) can 
change the signs of any four 5i for which the I\ lie in 
the same plane. Such transformations do not change the 
invariants (y\Viv$). By a sequence of these transforma- 
tions it is always possible to find a gauge in which Si = —I 
for at most one nonzero IV Define T* = Ti if there is 
one such point. If there is not, then set T* = 0. In this 
gauge, the mod 2 reciprocal lattice vector is G v = 2T*. 
The remaining invariant v$ is then determined by Si at 
Ti = 0. 

As we will explain below in section II. C. 2, the latter 
invariants, Vk are not robust in the presence of disor- 
der. We refer to them as "weak" topological invariants. 
On the other hand, vq is more fundamental, and distin- 
guishes the "strong" topological insulator. 

The formulas l2. 10112.121 are a bit deceptive because they 
appear to depend solely on a local property of the wave- 
functions. Knowledge of the global structure of |u ra k), 
however, is necessary to construct the continuous gauge 
required to evaluate Eq. 12.41 The existence of globally 
continuous wavefunctions is mathematically guaranteed 
because the Chern number for the occupied bands van- 
ishes due to time reversal symmetry. However, determin- 
ing a continuous gauge is not always simple. 



C. Surface States 

The spectrum of surface (or edge) states as a function 
of momentum parallel to the surface (or edge) is equiva- 
lent to the spectrum of discrete end states of the cylinder 
as a function of flux. Fig. 2 schematically shows two 
possible end state spectra as a function of momentum 
(or equivalently flux) along a path connecting the sur- 
face time reversal invariant momenta A Q and A 6 . Only 
end states localized at one of the ends of the cylinder 
is shown. The shaded region gives the bulk continuum 
states. Time reversal symmetry requires the end states 
at A a and A;, be twofold degenerate. However, there are 
two possible ways these degenerate states can connect up 
with each other. In Fig. 2a the Kramers pairs "switch 
partners" between A Q and Af,, while in Fig. 2b they do 
not. 

These two situations are distinguished by the Zi invari- 
ant characterizing the change in the time reversal polar- 
ization of the cylinder when the flux is changed between 
the values corresponding to A a and A 6 . Suppose that 
at the flux corresponding to A a the groundstate is non 
degenerate, and all levels up to and including the dou- 
blet e a i are occupied. If the flux is adiabatically changed 
to Af, then for Fig 2a the doublet Ebi is half filled, and 
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FIG. 2: Schematic representations of the surface energy levels 
of a crystal in either two or three dimensions as a function of 
surface crystal momentum on a path connecting A a and A5. 
The shaded region shows the bulk continuum states, and the 
lines show discrete surface (or edge) bands localized near one 
of the surfaces. The Kramers degenerate surface states at A a 
and At can be connected to each other in two possible ways, 
shown in (a) and (b), which reflect the change in time re- 
versal polarization ir a iTb of the cylinder between those points. 
Case (a) occurs in topological insulators, and guarantees the 
surface bands cross any Fermi energy inside the bulk gap. 
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FIG. 3: Diagrams depicting four different phases indexed by 
1/0; (U1V2V3). The top panel depicts the signs of Si at the points 
Ti on the vertices of a cube. The bottom panel characterizes 
the band structure of a 001 surface for each phase. The solid 
and open circles depict the time reversal polarization ir a at 
the surface momenta A a , which are projections of pairs of 
Ti which differ only in their z component. The thick lines 
indicate possible Fermi arcs which enclose specific A a . 



the groundstate has a twofold Kramers degeneracy asso- 
ciated with the end. For Fig. 2b, on the other hand, the 
groundstate remains non degenerate. This construction 
establishes the connection between the surface states and 
the bulk topological invariants. When 7r a 7r;, = — 1(+1) 
the surface spectrum is like Fig. 2a (2b). 

It follows that when 7r a 7Tj, = — 1(+1) a generic Fermi 
energy inside the bulk gap will intersect an odd (even) 
number of surface bands between A a and A b . Thus, when 
TTaTft = — 1 the existence of surface states is topologically 
protected. The details of the surface state spectrum will 
depend on the Hamiltonian in the vicinity of the surface. 
In Fig. 2 we have assumed that surface bound states 
exist for all momenta. This need not be the case, since 
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it is possible that by varying the surface Hamiltonian 
the degenerate states at A a and A b can be pulled out of 
the gap into the bulk continuum states. This, however, 
does not change our conclusions regarding the number 
of Fermi energy crossings. When 7r a 7r(, = —1 there still 
must exist surface band traversing the energy gap. 

In the two dimensional quantum spin Hall phase 
7Ti7T2 = — 1, and there will be an odd number of pairs 
of Fermi points 9 '^. In the simplest case where there is 
a single pair, the states at the Fermi energy will be spin 
filtered in the sense that the expectation value of the 
spin in the right and left moving states will have oppo- 
site sign. These states are robust in the presence of weak 
disorder and interactions because time reversal symme- 
try forbids clastic backscattering. Strong interactions, 
however, can lead to an electronic instability that opens 
a gap 31 ' 32 . The resulting ground state, however, breaks 
time reversal symmetry. 

In three dimensions, the Kramers degenerate band 
crossings that occur at A a in the surface spectrum are two 
dimensional Dirac points. While such Dirac points will 
occur in any time reversal invariant system with spin or- 
bit interactions, the nontrivial structure here arises from 
the way in which the Dirac points at different A a are con- 
nected to each other. This is determined by the relative 
signs of the four 7r a associated with any surface. 

In Fig. 3 we depict four different topological classes 
for three dimensional bandstructures labeled according to 
vq; {v\Viv%), which are characterized by different values 
of Si for the eight Ti , represented as the vertices of a cube 
in momentum space. The lower panel shows the surface 
Brillouin zone for a 001 surface with the four A a labeled 
by either filled or solid circles, depending on the value 
of ir a = 6i=(ai)5i=( a 2). The surface band structure will 
resemble Fig. 2(a) on paths connecting two filled circles 
or two empty circles, and will resemble Fig. 2(b) on 
paths connecting a filled circle to an empty circle. This 
allows us to draw a minimal Fermi arc (denoted by the 
thick lines) , which separates regions containing the filled 
circles from regions containing the empty circles. 



1. Strong Topological Insulator 

For the strong topological insulator, w± tt2^3^4 = — 1, 
so that a single ir a=a * differs in sign from the other three. 
The simplest Fermi arc, shown in Fig 3d thus encloses 
A Q * once. As in the two dimensional case, this Fermi arc 
can not be eliminated. In general, time reversal symme- 
try requires that any time reversal invariant Fermi arc C 
satisfy C — —C. It follows that the Berry's phase ac- 
quired by an electron circling the Fermi arc can only be 
either or n. Since the Fermi arc of the strong topologi- 
cal insulator encloses a single Dirac point an electron cir- 
cling the Fermi arc acquires a Berry's phase of ir. Since 
this can not be changed by continuous changes to the 
Hamiltonian, we conclude that the 7r Berry's phase is 
a generic feature of the surface Fermi arc in the strong 



topological insulator phase. Such a Fermi arc defines a 
"topological metal"—, which is topologically protected 
and, unlike an ordinary metal, can not be depleted. 

In the presence of weak disorder the ir Berry's phase 
changes the sign of the weak localization correction to 
the conductivity and gives rise to antilocalization, as in 
the symplectic universality clas o 16 i 18 . We suspect that in 
the absence of electron-electron interactions these surface 
states can not be localized even for strong disorder (pro- 
vided the bulk topological phase is not destroyed). As 
in the 2D case, however, electron interactions can lead 
to a localized phase, which will necessarily break time 
reversal symmetr y 22 ' 31 i 32 . 

In the strong topological insulator it is possible that 
the Fermi energy can be tuned to intersect a single Dirac 
point. This is a rather unique situation, because lattice 
Dirac fermions are generally expected to come in paira^ 9 -. 
These surface Dirac fermions are reminiscent of domain 
wall fermions which have been studied in the context of 
lattice gauge theories 3 - 3 -. The surface can be viewed as 
an interface between the topological insulator and a con- 
ventional insulator (the vacuum). These two phases can 
be characterized in terms of a three dimensional Dirac 
fcrmion, whose mass has opposite sign in the two phases 
(See for example section III.C). The domain wall between 
the two is then characterized by a gapless Fermion, or 
zero mode, which is related to the zero energy midgap 
states that appear in a one dimensional Pcicrls insula- 
tor at a soliton 34 . However, there are some important 
differences between our model and the conventional ap- 
plications of domain wall fermions. (1) In our problem 
there is no reason to have particle-hole symmetry, so tun- 
ing is required for the Fermi energy to be at the Dirac 
point. (2) The domain wall fermion applications have 
often been used to model chiral fermions in even dimen- 
sional space-time^ 3 -. Our 2+1 dimensional surface Dirac 
fermions are not chiral. Nonetheless, they realize the 2+1 
dimensional "parity anomaly"—. 

The parity anomaly arises for a single (i.e. undou- 
bted) species of massless Dirac fermion in 2+1 dimen- 
sions. When the response to the electromagnetic field is 
naively computed in this model, one finds 35 

J A = ±^ vX F IMV , (2.14) 

where J\ is the three current and is the electro- 
magnetic field tensor in 2 + 1 dimensions. This appears 
"anomalous" in the sense that the electromagnetic field 
gives rise to currents which appear to violate the sym- 
metries of the Dirac Hamiltonian. The sign ambiguity in 
(|2.14p is due to the regularization procedure, in which a 
finite mass is included to control divergences and taken 
to zero at the end. The origin of the singular behavior is 
the subtlety of this limiting procedure. 

In a magnetic field the Dirac equation leads to a Lan- 
dau level at exactly zero energy. At exactly half filling 
the system is thus at a critical point separating two quan- 
tum Hall states with a xy = ±(l/2)e 2 //i. This explains 
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the singular behavior described above. Indeed, the regu- 
lator mass term discussed above which opens a gap nec- 
essarily violates time reversal symmetry because it lifts a 
Kramers degeneracy. This leads to quantum Hall states 
even in zero applied magnetic field. 

For our problem, in the absence of time reversal sym- 
metry breaking perturbations we do not expect anoma- 
lous currents to occur. However, in a magnetic field, the 
parity anomaly shows up in the quantum Hall effect be- 
cause the surface Hall conductivity will be quantized in 
half integers, 

1 e 2 

a xy = (n+-) — . (2.15) 

It is interesting to ask whether such a "fractional" integer 
quantum Hall effect could be measured. Unfortunately, 
in a standard transport experiment in which currents and 
voltages are measured by attaching leads to a sample, the 
1/2 can not be directly measured. The problem is that 
in a slab geometry there is no simple way to isolate the 
"top" surface from the "bottom" surface. The two will 
always be measured in parallel, and the two half inte- 
gers will always add up to an integer. For such a trans- 
port experiment there is no getting around the "fermion 
doubling problem" . What is required is a method for 
measuring the current on the top surface without mea- 
suring the bottom surface. If it can be established that 
the currents are flowing on both surfaces, then dividing 
the measured Hall conductance by two could in principle 
demonstrate the half quantization. 

A lattice realization of the parity anomaly was pro- 
posed by Dagotto, Fradkin and Boyanovski 36,37 , w \ l0 
studied a tight binding model for PbTe in the presence 
of a domain wall where the Pb and Te atoms are in- 
terchanged. They showed that in their model the do- 
main wall exhibits massless Dirac fermions, and has a fi- 
nite Hall conductivity even at zero magnetic field. Their 
model, however is rather different from ours. In the pres- 
ence of the domain wall their Hamiltonian explicitly vi- 
olates time reversal symmetry^, and it leads to an even 
number of species of Dirac fermions on the domain wall. 
Haldane introduced a model of the quantum Hall effect 
on honeycomb lattice in a periodic mag netic field! This 
model, which also breaks time reversal symmetry realizes 
the parity anomaly (with a single Dirac fermion) when 
the Hamiltonian is tuned to the transition between the 
quantum Hall phase and the insulator. In this model, 
however, the Hall conductivity is an integer. 

The surface of the strong topological insulator is thus 
unique in that it can generate a single Dirac fermion with- 
out violating time reversal symmetry, and in principle 
exhibits the half quantized quantum Hall effect. 

2. Weak Topological Insulator 

When vq = 0, states are classified according to G„. 
We refer to the states with G v ^ as weak topo- 



logical insulators 15 , vq = implies that for any sur- 
face the associated time reversal polarizations will satisfy 
7ri7r27T37T4 = +1. This implies that either (1) all of the 
Tr' a s are the same or (2) two will be positive and two will 
be negative. The former case occurs for surfaces G = G„ 
mod 2, where G v is given in (|2.13|) . For these surfaces 
there are no topologically protected surface states. For 
G 7^ G„ mod 2, two of the A a 's are positive and two neg- 
ative. The Fermi arc encloses the two A a 's which have 
the same sign for 7r a . 

These states can be interpreted as layered two dimen- 
sional quantum spin Hall states. To see this, consider 
two dimensional planes in the quantum spin Hall state 
stacked in the z direction. When the coupling between 
the layers is zero, the electronic states will be indepen- 
dent of k z . It follows that the four Si's associated with 
the plane k z = n/a will have product —1 and will be the 
same as the four associated with the plane k z = 0. The 
topological invariants will then be given by vq = and 
G„ = (2,Tr/a)z. This structure will remain when weak 
coupling between the layers is introduced. More gener- 
ally, quantum spin Hall states stacked in the G direction 
will have G„ = G mod 2. This implies that quantum spin 
Hall states stacked along different directions Gi and G2 
are equivalent if Gi = G2 mod 2. 

The surface states can also be understood in this man- 
ner. When the coupling between the layers is zero, it is 
clear that the gap in the 2D system implies there will be 
no surface states on the top and bottom surfaces. On the 
sides, however, the Fermi points guaranteed for the edges 
of the 2 dimensional system will become straight Fermi 
lines, in the k z direction. The two Fermi lines will en- 
close two time reversal invariant momenta, which occur 
at k z = and k z — tt /a, as in Fig. [3) 

Since the surface Fermi arc encloses an even number 
of surface time reversal invariant momenta (and hence 
an even number of 2D Dirac points) it follows that the 
Berry's phase associated with the Fermi arc is zero. Thus, 
the surface states of the weak topological insulators do 
not enjoy the same level of topological protection as those 
of the strong topological insulator. Below we will argue 
that in the presence of disorder the weak topological in- 
variants lose their meaning. 



D. Z2 invariant and Boundary Condition 
Sensitivity 

Thouless, Niu, and Wu generalized the topological 
characterization of the integer quantum Hall effect to ex- 
press the Chern invariant in terms of the sensitivity of the 
groundstate of a bulk crystal to phase twisted boundary 
conditions 38 . This is more fundamental than the charac- 
terization in terms of Bloch wavefunctions because it does 
not rely on the translational symmetry of the crystal. It 
explains the topological stability of the Hall conductance 
in the presence of weak disorder. In this section we con- 
sider a corresponding generalization of the Zi invariant. 
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To do so, we consider large (but still finite) crystal with 
periodic boundary conditions in all but one direction. A 
phase twist e l9i is associated with each periodic boundary 
condition. This has the same structure as the cylinder 
(and generalized cylinder) considered in section IIB, but 
now the circumferences are much larger. The fluxes now 
correspond to the phase twists Oi = $ie/h. Since the 
cylinder is still finite the discrete states associated with 
the ends can be characterized by their degeneracy. This 
allows us to characterize the change in time reversal po- 
larization when the phase twists are changed by %. For 
non-interacting electrons, the invariants characterizing a 
large cylinder can be deduced from the bandstructure in- 
variants by treating the entire sample to be a unit cell of 
an even larger crystal. It is therefore necessary to con- 
sider the effect of enlarging the unit cell on the topological 
invariants. 

The 2D invariant v is preserved when the unit cell is 
enlarged. This is easiest to see by considering the effect 
of doubling the unit cell on the surface spectra of Fig. 
2. The effect of doubling the unit cell will be to fold 
the momenta A a and A& back on top of each other. If 
after enlarging the unit cell we add a small random po- 
tential, which lowers the translational symmetry, then all 
"accidental" degeneracies introduced by the zone folding 
will be lifted, while the Kramers degeneracies at the time 
reversal invariant momenta will remain. It is then clear 
that the manner in which the degenerate states at A a are 
connected to each other is preserved when the bands are 
folded back. Since the invariant v is unchanged when the 
unit cell is enlarged, it may be used to characterized the 
bulk system. 

In three dimensions the strong topological invariant vq 
is also invariant under enlarging the unit cell. This can be 
seen by noting that if the surface Fermi arc is folded back, 
the number of time reversal invariant momenta enclosed 
is preserved modulo 2. The weak topological invariants 
z/fe, however, are not preserved by this procedure. It is 
possible that for a Fermi arc which encloses two time 
reversal invariant momenta the two points can be folded 
back on top of each other. When the two bands are 
coupled to each other, a gap can then open at the Fermi 
energy, so that the Fermi arc can disappear. 

We thus conclude that the weak topological invariants 
are only defined for the bandstructure of a perfect crystal, 
and rely on the lattice translational symmetry. For a 
clean system, they have important implications for the 
surface state spectrum, but the topological distinction is 
eliminated by disorder. The strong topological invariant 
i/o, however, is robust, and characterizes a bulk three 
dimensional phase. 

The fragility of the weak topological invariants can also 
be understood by considering stacked two dimensional 
quantum spin Hall states. If two identical quantum spin 
Hall states are coupled together, the resulting two di- 
mensional system will be a simple insulator, and will not 
have topologcally protected edge states. Thus a stack of 
an even number of layers will be equivalent to a stack 



of insulators, while a stack of an odd number of layers 
will define a (thicker) 2D quantum spin Hall state. This 
sensitivity to the parity of the number of layers indicates 
that the weak topological invariants do not characterize 
a robust three dimensional phase, but rather, are only 
properties of the bandstructure. 

This formulation of the Z 2 invariant in terms of the 
change in the time reversal polarization of a large system 
as a function of twisted boundary conditions can be gen- 
eralized to interacting systems. This suggests that the 
strong topological insulator phase remains robust in the 
presence of weak interactions^. 

III. Z 2 INVARIANTS WITH INVERSION 
SYMMETRY 

In this section we show how the presence of inversion 
symmetry greatly simplifies the evaluation of the Z 2 in- 
variants. We will prove Eq. 1.1, which allows topological 
insulators to be identified in a straightforward manner. 

Suppose that the Hamiltonian Ti. has an inversion cen- 
ter at r = 0. It follows that [H, P] ~ 0, or equivalently 
H(— k) = PiJ(k)P _1 , where the parity operator is de- 
fined by 

P\r,s z ) = | -r,s z ). (3.1) 

Here r is the (3 dimensional) coordinate and s z is the spin 
which is unchanged by P because spin is a pseudovector. 
The key simplification for problems with combined inver- 
sion and time reversal symmetry is that the Berry cur- 
vature J- (It) must vanish because it must simultaneously 
be odd under time reversal k) = — .F(k)) and even 
under inversion k) = +JF(k)) 17 . The Berry curva- 
ture is given by Vk x A(h), where the Berry's potential 
is 

27V 

A{k) = -t 5^(Un )k |Vk|t»n,k}, (3.2) 

n=l 

where the sum is over the 2N occupied bands. The fact 
that = suggests it is possible to choose a globally 

continuous "transverse" gauge in which .4(k) = 0. We 
will show that in this special gauge the 5i defined in (|2.4p 
are given by (1-1), so the Z 2 invariants can be easily 
evaluated. 

The transverse gauge may be explicitly constructed as 
follows. In an arbitrary gauge consider the 2A^ x 2N 
matrix 

v mn (k) = (n m ,k|-P0|un,k)- (3.3) 

Since (a\b) = (86|Ga) and 9 2 = -1 it follows that u(k) 
is antisymmetric. Since [P<d,H(k)] = u(k) is unitary. 
Thus the Pfaffian of v(k) is defined, and has unit magni- 
tude. The phase of Pf[u(k)] depends on the gauge, and 
its gradient is related to .4(k): 

AQl) = ~TrKk)tV kW (k)] = -iV k logPf[u(k)]. (3.4) 
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The first equality is derived by differentiating (|3 . 3[) and 
using the unitarity of u(k). The second equality follows 
from the facts that det[u] = Pf[v] 2 and Vfelogdetfu] = 
Tr[V fc log v(k)} = Tr[«t(k)V k «(k)]. 

To set -4(k) = we thus adjust the phase of |w„k) to 
make 

Pf[u(k)] = 1. (3.5) 

This can be done, for instance, by a transformation of 
the form (j!T7]) . under which Pf[v(k)] -> Pf [u(k)]e" iek . 
In this gauge the problem of continuing \J det[w(k)] be- 
tween the Ti in ()2.4p is eliminated because det[u;(k)] = 1 
for all k. This can be seen by noting that v(k) has the 
property v(— k) = w(k)v(k)*w(k) T and using the iden- 
tity Pf[XAX T ] = Pf^jdetpf]. 

All that remains for Eq. (12. 4|) is to evaluate Pf[u>(Pi)]. 
To this end, we note that 

(ri) = (^ m , ri |P(Pe)|V„r i ). (3.6) 

Here we have used P 2 = 1, along with the anti- linearity of 
G, which allows us to replace |u n ri) by \ip n ri) = IV'n-r;) 
in Eq. (|2.5p . Since [H, P] = 0, \ip n ri) is an eigenstate of 
P with eigenvalue £, n (Ti) = ±1. After changing \ipnri) 
back to \unYi) it follows that 

(r,) = £ m (Ti)v mn (Ti). (3.7) 

The Pfafhan can be deduced from the following argu- 
ment, which uses the fact that the Pf [w] will be a poly- 
nomial in £„. First, note that 

Pf H 2 = det[iu] = det[v] Y[ Cn- (3.8) 

71=1 

Due the Kramers degeneracy, the distinct states \u2m.Ti) 
and l^m+i.r;} = ®\ u 2m,Ti) share the same parity eigen- 
value. Thus, each eigenvalue appears twice in the prod- 
uct in (|3.8[) . Taking the square root, we find 

N 

pf m = pf h n e 2m (3.9) 

The sign of the square root is fixed by the special case 
in which all £ n = 1, so that w = v. Since Pf[u] = 1 we 
conclude that in the transverse gauge, 

N 

$i=n 6m (r,). (3.io) 

m— 1 

Eq. (|3.10| is a central result of this paper. It means 
that with inversion symmetry the Z2 topological invari- 
ants can be deduced from the knowledge of the parity of 
each pair of Kramers degenerate occupied energy bands 
at the four (or eight in 3D) time reversal and parity in- 
variant points in the Brillouin zone. This provides a sim- 
ple method for determining the topological phase of any 



inversion symmetric insulator, without having to know 
about the global properties of the energy bands. 

In Eq. (|3.10p it appears as though each of the four 
(or eight) Si have gauge independent meaning, and thus 
provide extra topological invariants in addition to the 
one (or four) time reversal symmetry based invariants 
discussed in Section IIB. These extra invariants, however 
rely on the presence of inversion symmetry, and lose their 
meaning in the presence of surfaces, disorder or other 
perturbations which violate inversion symmetry. In con- 
trast, the invariants obtained from the product of 4 S^s 
do not rely on inversion symmetry for their existence. 
They depend only on time reversal symmetry, so they 
retain their value in the presence of inversion symmetry 
breaking perturbations. 



IV. TIGHT-BINDING MODELS 

In this section, we construct a class of inversion sym- 
metric tight-binding models that exhibit topological in- 
sulating states and apply the method presented in Sec- 
tion III to determine their topological classes. We will 
consider minimal models with 4 bands which result from 
four degrees of freedom per unit cell. We will focus on 
lattices in which the unit cell can be chosen to be inver- 
sion symmetric. We will see that this latter assumption 
makes the analysis of the topological phases particularly 
simple. While this assumption can always be satisfied 
for continuum models, it rules out certain inversion sym- 
metric lattice models, such as the rocksalt lattice. It is 
satisfied, however, for the specific examples we will con- 
sider. 

In IVA we study the general structure of this class of 
models, and then in IVB and IVC consider the specific ex- 
amples of the honeycomb lattice of graphene and the dia- 
mond lattice. In IVD we analyze a model for HgTe/CdTe 
quantum wells introduced recently by Bernevig, Hughes 
and Zhang^S. 



A. General Model 

We assume each unit cell associated with Bravais lat- 
tice vector R has four states |R, n). If the unit cell is 
parity invariant, then the parity operator P may be rep- 
resented as a 4 x 4 matrix as 

P\K,n) =^P nm | -R,m). (4.1) 

m 

In sections IVB and IVC we will consider examples in 
which each unit cell consists of two sublattices (denoted 
by Pauli matrix a z ) which are interchanged by inversion 
and two spin degrees of freedom (denoted by s z ). There- 
fore 

P = a x ®I, (4.2) 
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where I is the identity for the spin indices. In IVD P will 
have a slightly different form. The time reversal operator 
acting on the four component basis states is represented 
by 



9 = i(I <g> s v )K, 



(4.3) 



where K is complex conjugation, and / acts on the orbital 
indices. 

Given a lattice Hamiltonian TL in the |R, n) basis, we 
now consider the Bloch Hamiltonian 



H(k) 



ik-R 



(4.4) 



which for lattice periodic Bloch functions now becomes a 
4x4 matrix. Note that this transformation is slightly dif- 
ferent than the standard transformation of a tight bind- 
ing model for which R in (|4.4[) is replaced by r = R + d„ 
where d„ is a basis vector. The difference is a choice 
of gauge. With this choice H(k) has the properties 
H(k + G) = H(k) and H(-k) = PHi^p- 1 . Thus at 
the time reversal invariant momenta [H (k = 1^), P] = 0. 

It is convenient to express the 4x4 matrix H(k) in 
terms of the identity /, five Dirac matrices T a and their 
10 commutators r ab = [r a , T b ]/(2i) 41 . The Dirac matri- 
ces satisfy the Clifford algebra, r<T b + r<T a = 25 ab I. In 
this section, in order to avoid confusion of notation, the 
Dirac matrices r a will always appear with a superscript, 
and the time reversal invariant momenta will always be 
written as k = IV 

The choice of Dirac matrices is not unique. For ex- 
ample, in Ref. [ToL the Dirac matrices were chosen to be 
even under time-reversal, 0r°0 _1 — T a . In the presence 
of both inversion and time reversal symmetry it is more 
convenient to choose the Dirac matrices to be even under 
P0. Given the form of P and 0, the five matrices are 

r (i,2,3,4,5) = ^j^v® j ; a z g, s ^ fj z g, s y ; a * g s *), 

(4.5) 

With this choice of Dirac matrices the commutators are 
odd under P0, (P©^^©)- 1 = -r ab . Note that 
r 1 = P. It follows that 

o-pao-i _ ppap-i = 1 + r ° for a = 1 , . 

Ui U ri r \ -T a for a ^ 1. 

Time reversal and inversion symmetry imply that 
[P(k),P0] = 0. The most general Hamiltonian matrix 
is then 



H(k) =d a (k)I + Y,d a (.k)T a . 



(4.7) 



Written in this form, the energy eigenvalues naturally 
come in Kramers degenerate pairs with energy 



B(k) = d (k)±./X) d »( k ) 2 



(4.8) 




M 



M 



M 



n/2 
(b) 



FIG. 4: (a) Honeycomb lattice of graphene, with a unit cell 
indicated by the dashed lines, (b) Brillouin zone, with the 
values of Si associated with the time reversal invariant mo- 
menta labeled. Twa describes the loop enclosing half the zone 
used in Eq. l4~T3l 



At the time reversal invariant points k = T^, only T 1 
P is even under P and 0. Therefore 



H(k 



d (k = T i )I + d 1 (k = T i )P. 



(4.9) 



The parity eigenvalues £ ra for the states at k = Ti are 
given by the eigenvalues of P. It then follows from Eq. 
1.11 that provided there is an energy gap throughout the 
Brillouin zone, the Z2 invariants characterizing the va- 
lence band are determined by 



-sgn(d 1 (k = r i )). 



(4.10) 



We will use the above equation to determine the topolog- 
ical class of specific tight-binding models in the following. 



B. Graphene 

Graphene consists of a honeycomb lattice of carbon 
atoms with two sublattices. A tight-binding model which 
incorporates the symmetry allowed spin orbit interac- 
tions was introduced in Refs. andHo. 



H 



(ij) 



iX 



Gi 7 ' Cj 



(4.11) 



The first term is a nearest neighbor hopping term, and 
the second term is spin orbit interaction which involves 
spin dependent second neighbor hopping. This term is 
written in a way which can easily be generalized to three 
dimensions, s is the spin, and we have defined the unit 
vector 



4 X 4 



x d 2 .| 



(4.12) 



where dj- and df ■ are bond vectors along the two bonds 
the electron traverses going from site j to i. Thus, e«j -s = 
±s z . 

Choosing the unit cell shown in Fig. [4] the Hamilto- 
nian matrix H (k) can be determined using Eq. 14.41 and 
expressed in terms of Dirac matrices as in 14.71 The coef- 
ficients d a (k) are displayed in Table fl] The time reversal 
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do 







do 







di 


t [1 + COS X\ + COS X2] 




di 


t + Sti + t [COS X\ + COS X2 + COS X 3 ] 




d-2 


t[sina;i + sin £2] 




d 2 


t [sin xi + sin £2 + sin 353] 




d 3 







d 3 


Aso[sina;2 — sin 13 — sin(x2 — xi) + sin(a:3 


-xi)] 


di 







di 


\so[sinx 3 — sinii — sin(x3 — X2) + sin(xi 




d 5 


2Aso[sinxi — sinx2 — sin(a;i — X2)] 




d 5 


Aso [sin xi — sin x% — sin(xi — x 3 ) + sin(a;2 





TABLE I: Parameters for tight-binding model of graphene TABLE II: Parameters for diamond lattice tight-binding 
with xi = k ■ a; in a notation slightly different from Ref. Hol . model, with xu — k ■ a^ . 



invariant momenta, in the notation of Eq. 12.31 occur at 
k = r i=( - ni „ 2 ), for ri; = 0, 1. The Hamiltonian at these 
points can be deduced by noting that at k = Y i= ^ nin2 ^ 
xi = k • sli = riiir. r i= (oo) i s commonly referred to as 
the r point. The other three, which are equivalent un- 
der threefold rotations are called the M points. Using 
Eqs. (|1.1I1.2I4.10|) it is then straightforward to see that 
<5,:=(oo) = #i=(io) = *i=(oi) = -1j while 5 i= ( n ) = +1. The 
product is negative, so v = 1, and provided the energy 
gap is finite throughout the Brillouin zone the system is 
a topological insulator in the quantum spin Hall phase. 
The finite gap follows from the fact that e?o(k) — and 
there are no values of k for which all d a (k) = 0. 

The reason the three M points are not all the same is 
that the center of inversion defined by our unit cell is at 
the center of a bond, which does not have the threefold 
rotational symmetry. By choosing a different unit cell, 
with inversion center at the center of a hexagon, the M 
points would be equivalent. Our conclusion about the 
topological class, however, remains the same. 

It is interesting to note that the value of Si does not 
appear to have anything to do with the spin orbit in- 
teraction. The role that the spin orbit interaction plays 
is simply to ensure that the energy gap is finite every- 
where in the Brillouin zone. We will now argue for a 
parity and time reversal invariant system that if the spin 
orbit interaction is absent, then the negative product of 
Si implies that the energy gap must vanish somewhere in 
the Brillouin zone. This gives insight into the topological 
stability of the Dirac points in graphene in the absence 
of spin orbit interactions. 

We prove this by contradiction. In the absence of the 
spin orbit interaction we can consider spinless fermions. 
Suppose there is a finite gap everywhere, and the valence 
band is well defined throughout the Brillouin zone. Then 
on the one hand, the Berry curvature T = V x A is iden- 
tically zero due to inversion and time-reversal symmetry. 
On the other hand, we will show that the Berry's phase 
for the path t\i% shown in Fig. 0] which encloses half the 
Brillioun zone satisfies 

i£ A(k)-dk 

e ^1/2 = S 1 S 2 S 5 8 i . (4.13) 

Thus if 81828364 = —1, it would violate Stoke's theorem 
and leads to a contradiction. The tt Berry's phase thus 
requires that there either be a Dirac point in each half 
of the Brillouin zone, or a Fermi arc enclosing a Dirac 



point. 

To obtain Eq. (|4.13p for spinless electrons, we consider 
the unitary matrix 

my(k) = (iti,-k|P|%,k), (4.14) 

which is related to the Berry's potential via 
V fc logdet[m(k)] = -i(A(k) + A{~k)). Eq. (l4~T3D 
is then obtained by breaking the line integral into seg- 
ments connecting the time reversal invariant momenta 
and using the fact that det[m(k = I 1 *)] = Si. 

C. Diamond Lattice 

We now consider the tight binding model on a diamond 
lattice introduced in Ref. [lfj. This model exhibits both 
weak and strong topological insulator phases. 

The diamond structure consists of two interpenetrating 
face-centered cubic lattices displaced from each other by 
a basis vector d = a(l, 1, l)/4. The primitive translation 
vectors ai, a.2, a 3 are a(0, 1, l)/2, a(l, 0, l)/2, a(l, 1, 0)/2. 
Our model has the same form as Eq. 14.111 and includes a 
nearest neighbor hopping term as well as a second neigh- 
bor spin orbit interaction. 

It turns out that with this spin-orbit interaction term 
the valence bands and conduction bands meet at 3D 
Dirac points at the three inequivalent X points on the 
100, 010 and 001 faces of the Brillouin zone. In order to 
lift the degeneracy and obtain a gapped phase, we intro- 
duced a distortion, which changes the nearest neighbor 
hopping amplitudes. For simplicity we will focus here 
on a distortion in the 111 direction, which changes the 
nearest neighbor bond in the 111 direction, but leaves 
the other three bonds alone. The resulting model can be 
expressed in the form of Eq. 14.71 and the resulting d a (k) 
are listed in Table HT1 For Xso,8t ^ the gap is finite 
throughout the Brillouin zone. 

As in the previous section, the time reversal invari- 
ant momenta occur at k = r i= ( ni „ 2 „ 3 ) as in (|2.3[) . 
At these points xi = k • a; = n;7r. At the T point 
k = 0, {nin2n^) = (000). The three inequivalent X 
points (at k = (27r/a)(l, 0, 0) and related points) have 
(7x1712713) = (011), (101) and (110). The four inequiva- 
lent L points (at k = (7r/a)(l, 1, 1) and related points) 
have {n 1 n 2 n 3 ) = (100), (010), (001) and (111). The 111 
distortion makes the first three L points distinct from the 
fourth, which will be referred to as T. 
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From Table |TT] we can deduce the sign of d\ (k) , and 
hence Si at these points. We find Si = — 1 at T and the 
three L points, while <5, = +1 at T. At the X points 
5i = — sgn(Jii). Combining these we then find that 



(1; 111) for Stt > 
(0; 111) for Stt < 0. 



(4.15) 



When the 111 distorted bond is stronger than the other 
three bonds, so that the system is dimerized, the system 
is a strong topological insulator. When the 111 bond 
is weaker than the other three, so that the system is 
layered, it is a weak topological insulator with G„ = 
(27r/a)(l, 1, 1), which can be viewed as two dimensional 
quantum spin Hall states stacked in the 111 direction. 

In Ref. [l5j we computed the two dimensional band 
structure for the diamond lattice model in a slab geome- 
try. The results displayed the expected surface states, 
which behave according to the general principles dis- 
cussed in section II. C. 



D. Bernevig Hughes Zhang Model 

After this manuscript was originally submitted an in- 
teresting proposal appeared for the 2D quantum spin Hall 
effect in quantum well structures in which a layer of HgTe 
is sandwiched between crystals of CdTeML. Bernevig 
Hughes and Zhang (BHZ) showed that for an appropri- 
ate range of well thickness, the HgTe layer exhibits an 
inverted band structure, where the s and p levels at the 
conduction and valence band edges are interchanged. In 
this inverted regime, the structure exhibits a 2D quan- 
tum spin Hall effect. BHZ introduced a simple four band 
tight binding model which captures this effect. Though 
real HgTe does not have inversion symmetry, their toy 
model does. In this section we analyze this model and 
directly evaluate the Z 2 topological invariant using (1.1). 

BHZ considered a four band model on a square lattice 
in which each site has two Si/ 2 states \s, |), \s, J.) and two 
of the crystal field split p 3 / 2 states (with mj — ±3/2), 
\Px + ipyi T) an d \Px — ipy, [)■ The Hamiltonian is 

H = ^ £ a c\ aa C la<y - Y ta<T,ap4+a.<xa C iP° ( 4 - 16 ) 
i,<T,ot iaaaft 

where a labels the 4 nearest neighbors, er = ±1 describes 
the spin and a, (3 = s, p is the orbital index. The hopping 
term involves the matrix 



t. 



sp l - 



t p 

-t 



i<yO a 



(4.17) 



i>i> 



where 9 a gives the angle of nearest neighbor bond a with 
the x axis. 

As in Section IVA it is convenient to express this 
Hamiltonian in the form (|4.7p using Dirac matrices which 
are even under PQ. The form of the parity operator, 



do 


(e s + e P )/2 - 


(tss — t pp )(cosxi + COS £2) 


di 


(s s - s p )/2 - 


(t ss + t pp )(cosxi + cosa;2) 


d 2 




2t 3p sin x\ 


dz 







di 







ds 




2t 3p sina;2 



TABLE III: Parameters for the BHZ model, with xu = k • 



however is slightly different in this model, and Eq. 
replaced by 



P 



(4.18) 



where <r z = +1(— 1) describes s (p) states. The Dirac 
matrices are then the same as Eq. 14.51 except that a x 
and <r z are interchanged. The coefficients of these new 
Dirac matrices for this model are displayed in Table III. 

The analysis between Eqs. 14.61 and 14.101 proceeds ex- 
actly the same as before, and Si = — sgn[c?i(k = T,)]. We 
conclude that for k = (n/a)(ni, n 2 ), 



Si=(n 1 n 2 ) 



-sgn[ 



(t ss +u((-ir + M)" 2 )] 



2 v-sa ' u ppJ 

(4.19) 

For e s — e p > 4(< ss +t pp ) all of the (5 i= ( ni „ 2 ) are negative, 
so that the product v = +1. The system is a simple 
insulator. In this regime the bands have a conventional 
ordering throughout the Brillouin zone, with the s states 
in the conduction band and the p states in the valence 
band. For e s — e p < 4(t ss + t pp ) the bands near k = 
becomes inverted, and <5i = (oo) becomes positive, signaling 
a transition into the quantum spin Hall phase in which 
v = -1. 



V. TOPOLOGICAL PHASES IN SPECIFIC 
MATERIALS 

In this section we apply our method for evaluating the 
topological invariants to identify specific three dimen- 
sional materials that should exhibit a strong topological 
insulating phase. 



A. Bismuth Antimony Alloy 

Bi and Sb are group V semimetals in which there is 
a finite direct energy gap throughout the Brillouin zone, 
but a negative indirect gap due to band overlap. They 
have very close lattice parameters and form the solid al- 
loy Bii^Sb^i^. For .07 < x < .22 the indirect gap be- 
comes positive, leading to semiconducting behavior, with 
a maximum energy gap of order 30 meV for x — .18. In 
this section we will argue, based on the known band- 
structure of these materials, that this alloy is a strong 
topological insulator, which will have topological metal 
surface states. 



13 




.04 .07.09 .18 .22 

FIG. 5: Schematic representation of band energy evolution of 
Bii-zSbz as a function of x. Adapted from Ref. H^. 

Bulk bismuth and antimony have the rhombohedral A7 
structure, which consists of two interpenetrating, face- 
centered-cubic lattices which are displaced in the 111 di- 
rection and slightly distorted in the 111 direction. In 
bismut h 44 ' 45 , the valence band crosses the Fermi energy 
in the vicinity of the T point, which is located on the 111 
face of the Brillouin zone, giving rise to a small pocket 
of holes. The conduction band crosses the Fermi energy 
near the 3 equivalent L points, which reside at the other 
three body center zone faces, giving rise to pockets of 
electrons. At the L points, the bottom of the conduction 
band, which has L a symmetry is only slightly higher in 
energy than the next lower band, which has L a symme- 
try. In antimony^ 5 -, the electrons are again near the L 
point. However, unlike bismuth, the bottom of the con- 
duction band has L a symmetry. The holes are not at the 
T point, but rather at the lower symmetry H point. 

Despite the fact that bismuth and antimony have neg- 
ative indirect gaps, the finite direct gap throughout the 
Brillouin zone allows for the topological characterization 
of the valence energy bands. Since both bismuth and 
antimony have inversion symmetry, we can apply Eqs. 
(|l.lll.2|) by reading off the parity eigenvalues ^n(^i) from 
published band structure o 44 ' 45 . In Table IIVI we display 
the symmetry labels for the five occupied valence bands 
at the 8 time reversal invariant momenta (1^ = T, T, and 
the three equivalent L and X points). The parity eigen- 
value £ n (ri) can be read from the superscripts ± or the 
subscripts s/a = +/—■ (For an explanation of this nota- 
tion see Ref. lifil The right column displays the product 
Si from Eq. (|1.1[) . Based on this data, we conclude that 
the valence band of bismuth is equivalent to that of a con- 
ventional insulator, while the valence band of antimony 
is that of a strong topological insulator. The difference 
between the two is due to the inversion between the L s 
and L a bands. 

The evolution of the band structure of the alloy 
Bii^Sbj; as a function of x has been well studie d 42 ' 43 ' 47 , 
and is summarized in Fig. [5J As the Sb concentration 
is increased two things happen. First, the gap between 
the L s and L a bands decreases. At x = .04 the bands 



Bismuth 



ir 


r+ r~ r+ r+ r+ 

1 6 1 6 1 6 1 6 1 45 




3L 


L s L a L s L a La 




3X 


Xa X s X s X a X a 




IT 


rp — rp+ rp — rp+ rp — 

J 6 ± 6 J 6 J 6 J 45 






Z2 class 


(0;000) 


Antimony 


ir 


r+ r~ r+ r+ r+ 

1 6 1 6 1 6 1 6 1 45 




3L 


L s L a L s L a L s 


+ 


3X 


Xa X s X s X a X a 




IT 


T 6 - T+ T 6 ~ T+ T 4 - 5 






Z2 class 


(i; in) 



TABLE IV: Symmetry labels for the five valence bands of bis- 
muth and antimony at eight time reversal invariant momenta 
according to Ref. 1451 The parity eigenvalues can be read from 
+/— and s/a. Using Eqs (|l.lll.2|l they determine the topo- 
logical class. The indices (111) define a mod 2 vector (12.13[) 
in the direction of the T point. 



cross and the gap reopens with the inverted ordering. 
Secondly, the top of the valence band at T comes down 
in energy and crosses the bottom of the conduction band 
at x — .07. At this point the indirect gap becomes posi- 
tive, and the alloy is a semiconductor. At x = .09 the T 
valence band clears the L s valence band, and the alloy is 
a direct gap semiconductor at the L points. As x is in- 
creased further the gap increases until its maximum value 
of order 30 meV at x = .18. At that point the valence 
band at H crosses the L s valence band. For x > .22 the 
H band crosses the L a conduction band, and the alloy is 
again a semimetal. 

Since the inversion transition between the L s and L a 
bands occurs in the semimetal phase adjacent to pure 
bismuth, it is clear that the semiconducting Bii-^Sb^ 
alloy inherits its topological class from pure antimony, 
and is thus a strong topological insulator. Of course, this 
conclusion is predicated on a "virtual crystal approxima- 
tion" in which the disorder due to the random mixture is 
ignored, so that inversion symmetry is preserved in the 
alloy. However, since this inherent disorder does not de- 
stroy the bulk energy gap, it is unlikely to change the 
topological class, which does not require inversion (or 
translation) symmetry. We thus conclude that intrinsic 
Bii^Sbz, despite its bulk energy gap will have conduct- 
ing surface states, which form a topological metal. 

Semiconducting Bii-^Sba; alloys have been studied 
experimentally because of their thermoelectric proper- 
ties, which make them desirable for applications as 
thermocouple a 43 ' 48 ' 49 ' 50 . Transport studies have been 
carried out both on bulk samples^ 3 - and epitaxial thin 
filmss 9 -. For T > 50K semiconducting behavior is ob- 
served, while at lower temperatures the resistivity sat- 
urates at a value in the range 5 — 50/xf2m. This ob- 
served residual resistivity is probably too small to be ex- 
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FIG. 6: (a) Bandstructure of Q-Sn near the F point, which 
describes zero gap semiconductor due to the inverted Tg" and 
r^T bands, (b) In the presence of uniaxial strain, the degener- 
acy at F is lifted, opening a gap in the spectrum. The parity 
eigenvalues remain unchanged. 



Grey Tin 



ir 


1 6 1 7 


p- p+* 
1 7 1 8 








2Xs v 


+ 


4L 


A3 ^6 








z 2 


Class 


(1;000) 



TABLE V: Symmetry labels for the four valence bands of 
Grey Tin at eight time reversal invariant momenta according 
to Ref. Hi 

plained by surface states. It has been attributed to resid- 
ual charged impurities^, which act as shallow donors, 
making the alloy slightly n type. In order to separate 
the surface properties from the bulk transport, it will be 
necessary either to improve the purity of the samples, 
or perhaps use gating in a heterostructure to push the 
Fermi energy into the gap. 

B. Grey Tin and Mercury Telluride 

Tin is a group IV element, which in it's a (or grey) 
phase has the inversion symmetric diamond structure. 
Unlike carbon, silicon and germanium, though, it is a zero 
gap semiconductor, in which the ordering of the states 
at the conduction and valence band edge is inverte d 51 ' 52 . 
The Fermi energy lies in the middle of a four fold degener- 
ate set of states with symmetry, which can be derived 
from p states with total angular momentum j = 3/2. The 
four fold degeneracy at the T$ point is a consequence 
of the cubic symmetry of the diamond lattice. Apply- 
ing uniaxial strain lifts this degeneracy into a pair of 
Kramers doublets and introduces an energy gap into the 
spectrum 5 ^. For pressures of order 3 x 10 9 dyn/cm 2 , the 
induced energy gap is of order 40 meV. We now argue 
that this insulating phase is in fact a strong topological 
insulator. 

Table [V] shows the symmetry labels for unstrained 
a— Sn associated with the four occupied valence bands 
at the eight time reversal invariant momenta 5 -^. Uniaxial 
strain lowers the symmetry, so the cubic symmetry labels 
no longer apply. But since the strain does not violate in- 
version symmetry the parity eigenvalues are unchanged. 



The only effect is to split the degeneracy of the T% level 
into two sets of even parity Kramers doublets. In table 
IVl T^* refers to the occupied doublet. Based on the par- 
ity eigenvalues we conclude that strained grey tin is a 
strong topological insulator. 

HgTe is a II- VI material with the zincblend 
structur e 52 ' 55 . It is a zero gap semiconductor with an 
electronic structure closely related to grey tin. The Fermi 
energy is in the middle of the four fold degenerate Tg, 
states, whose degeneracy follows from the cubic symme- 
try of the zincblend lattice. As in grey tin, uniaxial strain 
lifts this degeneracy and opens a gap at the Fermi energy. 

Though HgTe lacks inversion symmetry, we now ar- 
gue based on adiabatic continuity that the gap induced 
by uniaxial strain leads to a strong topological insulator. 
The electronic structure of II- VI materials can be under- 
stood by adding an inversion symmetry breaking pertur- 
bation to a inversion symmetric group IV crysta l 52 ' 56 . 
Provided this perturbation does not lead to any level 
crossings at the Fermi energy, we can conclude that the 
II-IV material is in the same topological class as the 
group IV crystal. The bandstructures of grey tin and 
HgTe are very similar, and the cubic symmetry labels of 
the energy bands show how the bands evolve between the 
two. This allows us to conclude that strained a-Sn and 
HgTe will be in the same topological class, which is that 
of the strong topological insulator. 

In Ref. [251 Murakami. Nagaosa and Zhang introduced a 
four band tight binding model based onp 3 / 2 atomic levels 
on a fee lattice to describe strained a-Sn and HgTe. As 
argued in Ref. [l^, this model predicts that these mate- 
rials are simple insulators in the 0; (000) class. This can 
be understood by noting that since the model includes 
only p 3 /2 atomic levels the parity eigenvalues in (|1.1[) are 
all = —1. This contradicts the known band struc- 
ture of these materials, as displayed in Table [V] This 
model correctly describes the electronic states near the T 
point, but it gets the global topology of the bands wrong. 
To capture the global topology a tight binding model of 
these materials must include both s and p levels. The 
more recent theory^ of the 2D quantum spin Hall effect 
in HgTe/CdTe quantum wells discussed in Section IVD 
correctly incorporates s and p levels. 



C. Lead-Tin Telluride 

PbTe and SnTe are narrow gap IV- VI semiconductors 
with the rocksalt structure 5 ^. The bandgap in these ma- 
terials is direct, and occurs at the 4 equivalent L points 
in the FCC Brillouin zone. PbTe has an inverted band- 
structure relative to SnTe, in which the relative ordering 
of the Lq and L§ bands at the conduction and valence 
band edges are switched. Nonetheless, both of these ma- 
terials are conventional insulators. In Table IVIi we dis- 
play the symmetry labels at the 8 time reversal invari- 
ant points (r, 3 equivalent X points and 4 equivalent L 
points) 59 . Since the inversion occurs at an even number 
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TABLE VI: Symmetry labels for the five valence bands of Tin 
Telluride and Lead Telluride at eight time reversal invariant 
momenta, according to Ref. [Ssl 




FIG. 7: Schematic uniaxial strain-composition phase diagram 
for Pbi-zSn^Te. Away from the inversion transition at x ~ .4 
the material is a conventional insulator (I). Near the transition 
it is a strong topological insulator (STI). 



of points in the Brillouin zone, both materials belong to 
the conventional insulator topological class. 

PbTe and SnTe form an alloy Pbi-^Sn^Te. At x ~ .4 
there is an inversion transition where the band gap at 
the four L points vanishes, giving rise to three dimen- 
sional Dirac point o 57 ' 58 . The phases on either side of 
this transition are only distinct if inversion symmetry is 
present. Thus disorder, which is inevitably present in the 
alloy blurs the transition. However, uniaxial strain ori- 
ented along the 111 direction will distinguish one of the 
L points (call it T now) from the other three L points. 
It follows that the inversion transitions at the L and T 
points will occur at different values of x. Thus there will 
be an intermediate phase in which L is inverted, but T 
is not (or vice versa). From Eqs. (1.1,1.2) this interme- 
diate phase will be a strong topological insulator with 
surface states forming a topological metal. Note, that 
this direction depends on the orientation of the uniaxial 
strain. For instance strain along the 100 direction will 
distinguish two L points from the other two, and will not 
lead to an intermediate topological phase. In Fig. 2 we 
show a schematic phase diagram as a function of x and 
111 strain. 



The band inversion between SnTe and PbTe has been 
discussed by a number of authors previously. Volkov 
and Pankrato\*££ argued that PbTe and SnTe can be de- 
scribed by a low energy field theory consisting of three 
dimensional Dirac fermions with masses of opposite sign. 
They concluded that a planar interface between PbTe 
and SnTe will exhibit interface states described by a two 
dimensional massless Dirac equation. The appearance of 
such "domain wall fermions" is similar to the appearance 
of midgap states in a one dimensional Peierls insulator 
at a solitor*2i. A related proposal was made by Fradkin, 
Dagotto and Boyanovsk i 36 i 37 , who considered a domain 
wall in PbTe in which on one side the Pb and Te atoms 
are interchanged. This was also modeled as 3D Dirac 
fermions with a mass which changes sign at the inter- 
face. 

The domain wall fermions which appear in these the- 
ories are similar to the states we predict at the surface 
of a strong topological insulator. Indeed, if one views 
the vacuum as a band insulator with a large gap, then 
the surface can be viewed as an interface between a band 
insulator and a topological insulator, which can be de- 
scribed as an inversion transition, where there is a change 
in the sign of the mass of a 3D Dirac fermion. However, 
there is an important difference between the proposals 
discussed above and the surface states of the topological 
insulator: the strong topological insulator - band insu- 
lator interface involves a sign change in an odd number 
of Dirac points, while the interface models above involve 
four Dirac points. Having an odd number is crucial for 
the topological stability of the surface states. 



D. Other materials 

The materials we have proposed above should not be 
considered to be an exhaustive list. In general it is 
necessary to consider insulators composed of heavy el- 
ements. Another candidate for a topological insulating 
phase is Bi2Te3, which, like Bii-^Sb^, is known for its 
thermoelectric properties^. This material is also a nar- 
row gap semiconductor, with an energy gap of order .13 
eV. Though the crystal structure of this material is in- 
version symmetric, we have been unable to locate band 
theory calculations which display the parity eigenvalues. 

Another possible candidate is the zincblend semicon- 
ductor /3-HgS. The electronic structure of this material 
has been a subject of some controversy. According to 
Delhi 6 -, it is a semiconductor which has an unusual band 
ordering, with the and Tg levels in the valence band 
and the Tj level in the conduction band. If this is the 
case, we expect the material to be a strong topological 
insulator. However, this conclusion has been challenged 
by Moon, et al.— , who find a more conventional band 
ordering with the Tg level in the conduction band and 
the F7 and T% levels in the valence band, leading to a 
conventional insulator. 
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VI. EXPERIMENTAL IMPLICATIONS 

We now briefly consider possible experimental probes 
of topological insulators. We will focus here on the three 
dimensional strong topological insulator phase, for which 
we suggested several materials in the previous section. 

The most direct probe of the topological insulators is 
transport. Since there is a bulk excitation gap, trans- 
port in intrinsic samples at very low temperature will be 
dominated by the surfaces, which can be probed by the 
geometry dependence of the conductance. For example, 
for a wire geometry the conductance will be proportional 
to the circumference of the wire, rather than the area. 

In addition, since the topological metal phase is in 
the symplectic universality class the conductivity is ex- 
pected to increase logarithmically at low temperature: 
ct(T) cx \og{L in (T)/£)21, where L m is the inelastic scat- 
tering length and £ is the mean free path. 

An interesting prediction for the surface states is that 
due the tt Berry's phase associated with the surface Fermi 
arc, the surface quantum Hall effect in a perpendicular 
magnetic field should be quantized in half odd integers, 
o~ X y = (n+l/2)e 2 /h. As discussed in section II.C.l this is 
difficult to measure directly without separately measur- 
ing the currents flowing on the top and bottom surfaces 
of the sample. However if the parallel combination of the 
two surfaces could be measured, the resulting Hall effect 
would be quantized in odd multiples of e 2 /h. This is sim- 
ilar to the quantum Hall effect in graphen o 65 ' 66 , which 
is quantized in odd multiples of 2e 2 /h. The difference 
is due to the fact that graphcnc has four Dirac points, 
including spin. 

A practical difficulty with transport measurements is 
that it is necessary to distinguish surface currents from 
bulk currents. Since bulk currents scale with the sample 
width W, even though there is a bulk energy gap E g , 
the temperature must be low: T <C E g /\ogW/a, where a 
is the lattice constant. Moreover, since the materials we 
have suggested have rather small energy gaps, samples 
with high purity will be required to reach the intrinsic 
limit. As discussed in section IVa, the low temperature 
behavior of Bii^ Sb^ is dominated by a low concen- 
tration of charged impurities, which form an impurity 
band 43 . This is a ubiquitous problem for narrow gap 
semiconductors, due to their low effective mass and high 
dielectric constant. Clearly it would be desirable to have 
a transport geometry which probes the surface states, 
while controlling the position of the bulk Fermi energy. 
Perhaps this may be possible in a clever heterostructure 
geometry, where the bulk Fermi energy can be adjusted 
with a gate. 

An alternative probe of the surface states would be 
to map the surface Fermi arc using angle resolved photo 
emission. Such measurements could establish that the 
surface Fermi arc encloses an odd number of time rever- 
sal invariant momenta in the strong topological insulator 
phase. Detailed ARPES studies have been carried out 
on the surfaces of bismuth 67 ' 68 ' 69 and antimony 70 . How- 



ever, the presence of the bulk Fermi surface complicates 
the analysis of these materials. It would be interesting to 
see how the results are modified in the semiconducting 
Bii_a; Sba; alloy. 

Finally, since the surface states are spin filtered, electri- 
cal currents flowing on the surface will be associated with 
spin accumulation, leading to a spin Hall effect. In GaAs, 
spin accumulation on a surface has been measure d 71 ' 72 . 
The narrow energy gaps in our proposed materials make 
detection of the spin accumulation more difficult. Per- 
haps a heterostructure geometry could make this possi- 
ble. 



VII. CONCLUSION 

In this paper we discussed topological insulator phases 
in two and three dimensions. We discussed in detail 
how the Z2 topological invariants characterizing these 
phases influence the surface state spectrum. In particu- 
lar, the quantum spin Hall phase in two dimensions and 
the strong topological insulator in three dimensions have 
robust conducting surface states, which we have char- 
acterized as a topological metal. We showed that the 
Z2 invariants can be determined easily in parity invari- 
ant crystals from the knowledge of the parity eigenval- 
ues for states at the time reversal invariant points in the 
Brillouin zone. Using this method, we deduced that the 
semiconducting alloy Bii_. c Sh r is a strong topological in- 
sulator, as are a-Sn and HgTe in the presence of uniaxial 
strain. 

There remain a number of further issues which need to 
be understood better. High among them are the effects 
of disorder and interactions. These are important both 
for the topological metal surface states as well as for the 
bulk topological phases. Numerical work by Onoda et 
ali^ 3 . has suggested that the transition between the con- 
ventional insulator and the quantum spin Hall phase in 
two dimensions belongs to a new universality class. It 
will be of interest to understand this transition better, 
along with the related transition between the topological 
insulator and the Anderson insulator, which presumably 
occurs when disorder is increased beyond a critical value. 

Finally, it would be desirable to develop a field the- 
ory for the topological insulating phases analogous to the 
Chern Simons theory of the quantum Hall effect. Perhaps 
this may lead to analogs of the fractional quantum Hall 
effect for the topological insulators. 
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